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Abstract 

We propose a new scenario in a radiative seesaw model based on flavor symmetry. In this 
model, we explore a possibility of linking non-zero #13 to dark matter. And we analyze the lepton 
sector to predict the observed neutrinos and mixings, especially obtaining a lower bound of $13 > 
3.5°. We show that the non-zero #13 is correlated with our heavy Majorana type of dark matter. 
Also we predict that the mass be 0(1-10) TeV, as a result of analyzing the Wilkinson-Microwave- 
Anisotropy-Probe and lepton flavor violation. 
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I. INTRODUCTION 



The large values of the solar (#12) and atmospheric (#23) mixing angles in the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix [l| may be telling us about some new symmetries of 
leptons not presented in the quark sector and may provide a clue to the nature of the quark- 
lepton physics beyond the standard model (SM). If there exists such a flavor symmetry in 
Nature, the tribimaximal (TBM) j^j pattern, sin 2 6 , 12 = 1/3, sin 2 # 23 = 1/2, sin #13 = 0, for 
the neutrino mixingwill be a good zeroth order approximation to reality. A Non-Abelian 
discrete symmetry [3| could play an important role in predicting the mixing angles. For 
example, in a well-motivated extension of the SM through the inclusion of A4 discrete 
symmetry, the TBM pattern comes out in a natural way in the work of flfl. Although 
such a flavor symmetry is realized in Nature leading to exact TBM, in general there may be 
some deviations from TBM [7, 8|. Recent data of the T2K, MINOS and Double Chooz 9] 
Collaborations imply that the unknown mixing angle, #13, can be relatively large, indicating a 



deviation from the exact TBM which 
And, the analysis based on global fits 



eads to vanishing CP violation in neutrino oscillations. 



10Nl2| of neutrino oscillations enter into a new phase of 



precise measurements of the neutrino mixing angles and mass-squared differences, indicating 
that the TBM mixing for three flavors of leptons should be modified. Their current best-fit 
values in la (3a) ranges from neutrino oscillation experiments are given by [lo| 
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which corresponds to normal neutrino mass ordering 1 . And we know nothing about all 
three CP-violating phases Sep, S± and 62- Besides the mystery of the mixing pattern, tiny 
neutrino mass is one of the most challenging problem beyond SM. It suggests that neutrinos 
could be induced by a radiative correction or tiny couplings if a theory be within TeV scale. 
Several years ago, Ernest Ma introduced the so-called radiative seesaw mechanism 
the masses are generated through one- loop effects [jjj]. Moreover, since the existence of the 
flavor nentrino mixing for the three neutrinos ^ imp.ies that the indiv.dnal lepton 
charges, L a , (a = e, /i,r), are not conserved [15J, the observation of neutrino oscillation 
have the possibility of measurable branching ratio for charged lepton flavor violation (LFV) 



131 ] where 



See Ref. 



10 ] for the data of inverted mass ordering. 
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decays such as /x — > ery, r — > cy and r — > /17, etc.. Experimental discovery of such lepton 
rare decay processes is one of smoking gun signals of physics beyond SM. 

Our starting point is an effective Lagrangian obeying the discrete A 4 x Z 2 symmetry which 
is spontaneously broken by the vacuum expectation value (VEV) of SU(2)l x U(1)y singlet 
scalar fields at a scale higher than the electroweak scale. In addition, we assign a Z 2 -odd 
quantum number to a new Higgs doublet and three right-handed singlet fermions while all 
SM particles are Z 2 -even parity, in order to explain both the TBM at tree level Lagrangian 
and dark matter (DM). After electroweak symmetry breaking, the Z 2 symmetry is exactly 
conserved and a Z 2 -odd doublet Higgs that does not have VEV, so called inert Higgs, while 
the standard Higgs boson get a VEV, which means the Yukawa coupling corresponding to 
Z 2 -odd Higgs doublet will not generate the Dirac mass terms in neutrino sector. Thus, the 
usual seesaw mechanism does not work any more and we naturally have a good candidate 
of DM corresponding to the lightest Z 2 -odd particle or Large Hadron Collider (LHC) sig- 
nals through the standard gauge interactions in our scenario. The assigned leptonic flavor 
symmetry will lead us to a neutrino mass matrix through one-loop mediated by a new Higgs 
doublet and right-handed neutrinos having Z 2 -odd parity, and indicating the deviations from 
TBM in lepton sector by dimension-5 effective operators driven by SU(2)l x U(1)y singlet 
scalar fields with a cutoff scale A. Here we assume that CP is a good symmetry above the 
cut-off scale A in Lagrangian level. In this paper, we address the possibility of a linking 
between DM and a non-zero 6*13 through the combination of WMAP results 16] with LFV 
/jL — y cy decay. We analyze possible spectrums of light neutrinos and their flavor mixing 
angles. In order to investigate the relation between DM and a non-zero #13, we focus on a 
normal hierarchical mass spectrum of light neutrinos. 

This paper is organized as follows. In section II, we show our model building for the 
lepton sector, in which we adopt A4 group. In section III, we discuss the predictions for 
neutrinos coming from the flavor symmetry. In section IV, we show constraints of DM mass 
from WMAP and LFV; \i —> ej, r — > ej, and r — > /ry processes. Section V is devoted to 
conclusions/ discussions. We discuss the Higgs potential in the appendix. 
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II. FLAVOR A A SYMMETRY AND A DISCRETE SYMMETRY Z 2 



In the absence of flavor symmetries, particle masses and mixings are generally unde- 
termined in gauge theory. We work in the framework of the SM, extended to consist of 
the right-handed S"?7(2)£,-singlet Majorana neutrinos, Nr. The scalar sector, apart from the 
usual SM Higgs doublet $, is extended through the introduction of two types of scalar fields, 
X and ?7, that are singlet and doublet under SU{2) L x U(l)y, respectively: 

$=(^°) T , X , V={V + V°) T - (2) 

Furthermore, to understand the present neutrino oscillation data we consider flavor 
symmetry for leptons, and simultaneously for both the existence of DM and TBM at tree- 
level Lagrangian to be explained we also introduce an auxiliary discrete symmetry Z 2 in 



a radiative seesaw [13] . Here we recall that is the symmetry group of the tetrahedron 
and the finite groups of the even permutation of four objects. Its irreducible representations 
contain one triplet 3 and three singlets 1,1', 1" with the multiplication rules 3 © 3 = 3 S © 
3 a © 1 © T © 1", T © 1" = 1, T © T = 1" and 1" © 1" = T. Let us denote (ai, a 2 , a 3 ) and 
(b\, b 2 , 63) as two A4 triplets, then we have 

(a©6) 3s = (a 2 b 3 + a 3 b 2l a 3 bi + ai& 3 , a x b 2 + a 2 b l ) , 
(a © 6) 3a = (a 2 b 3 - a 3 b 2 , a 3 6i - aib 3 , a x b 2 - a 2 b x ) , 
(a © 6)1 = aibi + a 2 b 2 + a 3 b 3 , 



(a © b) v = ai&i + uja 2 b 2 + u a 3 b, 
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(a©6)i» = ai&i + u) a 2 b 2 + ooa 3 b 3 , (3) 

where a; = e l2?r / 3 is a complex cubic-root of unity. 

The field contents under SU(2) x U(l) x A 4 x Z 2 of the model are assigned in Tabled: 

TABLE I: Representations of the fields under A4 x Z 2 and SU(2)l x U(l)y- 
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At the Lagrangian level, we shall assume the absence of CP violation in the Dirac neutrino 
sector and in charged lepton Yukawa interactions above the cutoff scale A, which for scales 
below A is expressed in terms of effective dimension-5 operators. With dimension-5 operators 
driven by X fields the Yukawa interactions (d < 5) in the neutrino and charged lepton sectors 
invariant under SU(2) x U(l) x A4 X Z 2 can be written as 

-£ Yuk = y v (L L N R ) 1 ij + ^M(I^N R ) 1 + ^(MNR)a.X 

+ y^[(ii^ R )^fi+^[(L L N R )^fi 

+ y e (L L $)il R + y^L L §) v l' R + y T {L L <&) v ,l' R 

+ V f[{L L ^h sX ]l R + V ^[(L L ^ 3sX ]vl' R + y f[(L L <S>) 3sX ]r>l' R 

+ yl[(L L ^) 3aX ]l R + V l[(L L <f>) 3aX ] v l R + V f[{L L <S>) 3aX ] v ,l' R + h.c , (4) 

where fj = 2r 2 ?7* and L L = [y L ,t^) T are the Higgs doublet and the lepton doublet trans- 
forming as singlet 1 and triplet 3 under A±, respectively, and here r 2 is the Pauli matrix. In 
the above Lagrangian, heavy neutrinos acqiure a bare mass term M and a mass induced by 
the electroweak singlet X scalar with 3 representation under A±. By imposing an additional 
symmetry Z 2 as shown in Tabled the A 4 x SU(2)l X U(1)y invariant Yukawa term £lN r $> 
is forbidden from the Lagrangian, and the neutral component of scalar doublet rj will not 
generate a VEV, 

(V°) =v v = 0. (5) 

Therefore, the scalar field t] can only couple to the standard gauge bosons as well as the Dirac 
neutrino mass terms are vanished, which means the usual seesaw does not operate anymore. 
However, the light neutrino Majorana neutrio mass matrix can be generated radiatively 
through one- loop with the help of the Yukawa interaction L^Nnf] and LiN^-fj/A in the 
Lagrangian, we will discuss it more detaily in Sec. III. 

Taking the A4 symmetry breaking scale above the electroweak scale in our scenario, that 
is, (x) > ( < ^ )0 ), and assuming the vacuum alignment of fields (%») as 

(Xi) = v x ^ 0, ( X 2> = (X3> = , (6) 
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then the right-handed neutrino Majorana mass terms are given by 



M R = M 



( 



1 

1 ne^ 
ne* 1 



\ 



(7) 



where k = |A* v x /M\ with (xi) = v Xi (i = 1,2,3). After the electroweak symmetry breaking 
(77°) = v with the VEV alignment in Eq. ([6]), the Yukawa interaction y^LiNnfj, together 
with the terms ^-[(L L N R ) 3a x\ • fj can be written as vVlY u Nr with the neutrino Yukawa 



coupling matrix Y v given by 



e ip \y,\ 



1 



y 2 e 



Vie 



i in 



t/J-2 



(8) 



/ 



where yi = \y s N +y%\v x /\y u \A, y 2 = \y s N -y^\v x /\y u \A and p = aig(y u ). Eq. (ED indicates that, 
once the VEV alignment in Eq. OH]) is taken, the A 4 symmetry is spontaneously broken and 
its residual symmetry Z 2 js, 7] is also explicitly broken by the higher dimensional operators. 
Therefore, we can expect a low energy CP violation responsible for the neutrino oscillation as 
well as high energy CP violation responsible for the leptogenesis in the neutrino sector, which 
could be generated by the off-diagonal terms of Yukawa neutrino coupling matrix [3, 

Assuming equally aligned VEVs of triplets, ($°) = (v,v,v), together with the VEV 
alignment in Eq. fl^]), the light charged lepton mass matrix can be explicitly expressed as 
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(9) 



where U W V[ and indicate that the left and right-diagonalization matrices for the charged 
lepton sector, and the mass matrix elements are written as 



mil = v (Vh + 2 V3), m i2 = 2vh 2 /3, 



rn 



13 



2vh 3 /3 
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21 



v{gi-hi)/3, m e 22 = v(y h + (g 2 -h 2 )/3), m 23 = v (g 3 - h 3 )/3 



m 



31 



-v(gi + hi)/3, m 



42 



-v(g 2 + h 2 )/3, 



m e 33 = v(y h -(g 3 + h 3 )/3) ,(10) 
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with hi = v x y s f JA, h 2 = v x y s f JA, h 3 = v x y s f J A, g x = -iV3v x y a f J A, g 2 = -iy/3v x y a f J A, 
gs = — iy/3v x yj 3 / A. All of them are complex. The important point is that, taking the 
VEV alignment of (xi) in Eq. fl6]) together with the equal VEV alignment of ($°), the A4 
symmetry is broken and its residual C 3 symmetry [5] is also broken through the dimension-5 
operators. One of the most striking features in the charged lepton sector is that the mass 
spectra of charged leptons are strongly hierarchical, i.e. the third generation fermions are 
much heavier than the first and second generation fermions. From Eq. OH]), for the most 
natural case of hierarchical charged lepton Yukawa couplings y/ 3 3> y/ 2 3> y^, the corrected 
off-diagonal terms which are from dimension-5 operators are not larger than the diagonal 
ones in size. Then, V[ and can be obtained by diagonalizing the matrices U^m£m\U u 
and mime, respectively, in Eq. (jHJ). Especially, the mixing matrix V[ takes part in PMNS 
mixing matrix. The matrix V[ can be, in general, parametrized in terms of three mixing 
angles and six CP violating phases : 



I c 2 c 3 c 2 s 3 e^3 s 2 e^2 \ 

— cis 3 e _l *3 — s 1 s 2 c 3 e l ^i - '^ C1C3 — s 1 s 2 s 3 e i ^~ (( ' e2+ ^ sic 2 e^i 
\ Sis 3 e~ l ^ + ^ — cis 2 c 3 e~^2 — SlC ^ e - l <t>\ _ Cl S2S 3 e l( -^~^ C\C 2 j 



Pi , (11) 



where Sj = sin^, q = cos6*j and a diagonal phase matrix Pe = Diag.(e*^, e^ 2 , e 1 ^) 
which can be rotated away by redefinition of left-charged fermion fields. In the charged 
fermion (quarks and charged leptons) sector, there is a qualitative feature which distin- 
guishes the neutrino sector from the charged fermion one. The mass spectrum of the 
charged leptons exhibits a similar hierarchical pattern as that of the down-type quarks, 
unlike that of the up-type quarks which shows a much stronger hierarchical pattern. For 
example, in terms of the Cabbibo angle A = sin^c ~ \V US \, the fermion masses scale as 
(m e , m M ) ~ (A 5 , A 2 ) m T , (m^, m s ) ~ (A 4 , A 2 ) m& and (m u ,m c ) « (A 8 , A 4 ) rru. This may lead 
to two implications: (i) the Cabibbo-Kobayashi-Maskawa (CKM) matrix 18] is mainly gov- 
erned by the down-type quark mixing matrix, and (ii) the charged lepton mixing matrix is 
similar to that of the down-type quark one. Therefore, we shall assume that (i) Vqkm = Vl 
and V£ = I, where V£ (V£) is associated with the diagonalization of the down-type (up-type) 
quark mass matrix and I is a 3 x 3 unit matrix, and (ii) the charged lepton mixing matrix 
V[ has the similar structure as the CKM matrix. Now, putting a reasonable assumption 

m„ rrinJ \Tm-A mi, mS, \m\o\ , s 

l™3sl Ksl \ m M l m 33l l m 33l Kssl I "^33 1 



into Eq. (J9]), then the measured mass hierarchy of charged lepton is expressed as m^/m T m 
\ m 22/ m z3\ ~ A 2 ,m e /m T ~ l^ii/Ti^l ~ 0.6A 5 and m e /m M « l m n/ m 22l ~ 0.5A 3 . And the 
mixing angles and phases in Eq. fill I) can be roughly expressed as 
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Now letting |(V/)i2| = IC2S3I ~ A in similar to the case of quark sector [7|, we obtain 



L)V1\ 



\m 



12I 



22 1 



A 2 \m 



121 



33 1) 



leading to |^i2l/l m 33l ~ A 3 . Consequently, 9\ 



|m 23 |/|m3 3 | ~ A 2 and 9 2 ~ I m i3l/I m 33l — A 3 are obtained from Eq. (1121) . Then, the mixing 
matrix V[ in Eq. (TTTj) can be written as 
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where the coefficients A and .B are real and positive coefficients, but less than one in mag- 
nitude. 

The Yukawa interactions in Eq. (J3J) and the charged gauge interactions in a weak eigen- 
state basis can be written as 



C = ^N c R M R N R + £ L mdR + WY„N R r) + ^W-£ L ^VL + h.c 



(15) 



When dealing with lepton flavor violation and DM it is convenient to work in the basis where 
heavy Majorana neutriino and charged lepton mass matrices are diagonal. In order to go 
into the physical basis (mass basis) of the right-handed neutrino, performing basis rotations 

L L ^V^UlL L , in^V^R, N R —> U R N R (16) 

where L L = (vl,£l) T , so that the Yukawa coupling matrix Y u gets modified to 

Y V ^Y V = V?UXUr , (17) 

and both the right-handed Majorana mass matrix M R and the charged lepton matrix me 
become diagonal by the unitary matrix Ur and U^Vl, respectively; 

m e = V?UlmiV& = Diag.(m e , m M , m T ) , M R = U^MrUr = MDiag.(a, 1, b) (18) 
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FIG. 1: One-loop generation of light neutrino masses. 



where a = a/1 + k, 2 + 2k cos£ and 6 = a/1 + k 2 — 2k cos £ with real and positive mass 
eigenvalues, Mi = Ma, M 2 = M and M 3 = M6, and the diagonalizing matrix Ur is 

1 



with the phases 



u 



R 



1 -1 
1 1 



1 



V 



e*-* 



/ 



tan 1 



Ksin£ 
+ k cos £ 



and tp2 = tan" 



ft sin £ 
« cos £ — 



(19) 



(20) 



III. LOW ENERGY OBSERVABLES 



We now proceed to discuss the low energy neutrino observables. Due to the Z2 symmetry, 
we can not get the Yukawa Dirac neutrino mass matrix and therefore the usual seesaw 
does not operate any more, however, similar to 19| the light neutrino mass matrix can be 
generated through one loop diagram in Fig. [T] due to the quartic scalar interactions. After 
electroweak symmetry breaking, i.e. (<&°) = vl, in a basis where charged lepton mass matrix 
is real diagonal the flavor neutrino masses can be written as 
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with Zi = Mf/mi, if mu(mi) is the mass of rf^ijfj) and m 



R(I) 



v,2 



± where 

the subscripts R{I) indicate real (imaginary) component 2 , respectively. With Mr = 
Diag(M rl , M r2 , M r3 ) and M ri = M i f~ 1 (z i ), the above formula Eq. (|2"TI) can be expressed 
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where rrii{i = 1, 2, 3) are the light neutrino mass eigenvalues and 
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As can be seen in Eq. (123j) . the leptonic mixing matrix is given as 

t/pMNS = , 

with the additional mixing matrix U v which is given by 

(i o oW 
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where the phase matrix of Majorana particle P u = Diag.(e l< ' 1 , e JI > 2 , e^ 3 ) can be absorbed into 
the neutrino mass eigenstates fields. The phase <^ 12 and the mixing angle 9 are given by 

2\AG* + GB*\ 



ip l2 = i p 1 -ip 2 = arg {AG* + GB*} 



tan 29 



(27) 



\B\ 2 - \A\- 

indicating that, in the limit of yi t2 approaching to zero, the angle 9 goes to vr/4(— 7r/4) and the 
phase ipi 2 goes to 7r(0) for - ^ < - ^ ( j ^ > ^ jp j > whose one-to-one correspondence 
comes from the constraint of mixing parameters presented in Eq. ( l37|) with experimental 
data, especially see the parameter e in Eq. 



See Appendix. 
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The light neutrino mass eigenvalues are given as 

ml = m 2 (|Gf + \A\ 2 cos 2 9 + \B\ 2 sin 2 9 - \AG* + OB* | sin 20) , 
m l = rn 2 f 2 (z 2 ) , 

m 2 3 = m 2 (\G\ 2 + \A\ 2 sin 2 9 + \B\ 2 cos 2 9 + \AG* + GB*\ sin 29) . (28) 

Because of the observed hierarchy |Am^ tm | ^ Amg ol > 0, and the requirement of Mikheyev- 
Smirnov-Wolfenstein resonance for solar neutrinos, there are two possible neutrino mass 
spectrum: (i) mi < m 2 < m 3 (normal mass spectrum) which corresponds to 9 = tt/4 + 5 
and (fi2 = 71 + 5', and (ii) m 3 < mi < m 2 (inverted mass spectrum) which corresponds to 
9 = — 7r/4 + 5 and <pi 2 = where \6\, \8'\ 1. The solar and atmospheric mass-squared 
differences are given by 

Am 2 Sol = m 2 2 - ml = m 2 {f 2 (z 2 ) - \G\ 2 - \A\ 2 cos 2 9 - \B\ 2 sin 2 9 + \AG* + GB*\ sin29} , 
A< =ml-ml = 2ml'- — ^ [ , (29) 

which are constrained by the neutrino oscillation experimental results. Note here that the 
parameter M r i can be simplified in the following limit cases as 



Mi [In Zi-l] 1 , for z { > 1 

2M U for Zi -> 1 (30) 



m 2 M~ l , for ^ < 1 . 



In this work we will focus on the case the lightest Z 2 -odd neutral particle of Ni which 
is stable and can be a candidate of DM. As will be seen in section- HVt DM constraint 
(Br(/i — > e7)+WMAP analysis) gives m v ± pa Mi ig htest where Mi ig htest is the lightest of the 
heavy Majorana neutrinos. And, due to lA^I < 47r, one can obtain Miightest ~ m rj (see 
Appendix). From Eq. (fl8|) . depending on the hierarchy of the heavy Majorana neutrino 
masses M 1 ,M 2 and M 3 , the relative size of the parameter k consistent with the possible 
mass ordering of light neutrinos and hierarchy of Am^ and Am^ can be classified as follows 

(i) M 3 < M 2 < Mi for f = 0: 

• M 3 ~ m v M 2 < Mi (or a > 1 3> b) : this case corresponds to the normal 
hierarchical mass spectrum with b — > i.e. k ~ 1 and a ~ 2. Using f(zi) ~ 
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21n(2/6) - 1 > f{z 2 ) ~ -2 In 6 - 1 > f(z 3 ) ~ 1/2, the condition Am 2 21 > is 
obtained for = +7r/4, indicating the normal mass ordering of light neutrinos with 
Am 32 « m 2 , (f 2 (z 3 )/b 2 — f 2 (zi)/A) > 0. The ratio of the mass squared differences 
defined by R = Am^/IAmSy is given by 

i2 « 6 2 (4/ 2 (^ 2 ) - / 2 (^)) , (31) 

where the equality roughly is given under y 1>2 , b <C 1. From Eq. (j3T|) and the best-fit 
values of the solar and atmospheric mass squared differences (R ~ 3 x 10~ 2 ), one can 
determine the magnitude of the parameter b, roughly 

&» 0.01 . (32) 

• M 3 ~ fhr! < M 2 < Mi (or b < 1 < a) : this case gives a ~ 1.4, 6 ~ 0.6 with n ~ 0.4 and 
a ~ 2.8, b ~ 0.8 with k ~ 1.8 in numerical calculations, which corresponds to f(z 3 )/b < 
f{z\)/a < /(^) indicating a degenerate inverted orderings of light neutrinos. 

(ii) Mi ~ < M 2 < M 3 for £ = 7r : this case a < 1 < 6 corresponds to a degenerate 
inverted ordering of light neutrinos giving a ~ 1 and 6 ~ 3 with k ~ 2 in numerical 
calculations. Note that a case Mi ~ <C M 2 < M 3 is not allowed, due to e — > 2 in 
Eq. (1371) . actually because it could not satisfy the experimental data of light neutrino of 
mixing angles. 

(iii) M 2 ~rh v < M x < M 3 (1 < b < a) for £ = : this corresponds to f(z 1 )/a < f(z 2 ) < 
f(z 3 )/b, giving a ~ 4.7 and b ~ 2.7 with k ~ 3.7 in numerical calculations. It gives a 
degenerate normal mass ordering of light neutrinos. Note that a case M 2 ~ fh v <C Mi ~ M 3 
for £ ~ 7r/2 or £ ~ 37r/2 is not allowed because it could not satisfy the ratio R ~ 3 x 10 -3 . 

Note here that in our scenario the inverted hierarchical light neutrino mass spectrum is 
not allowed because the condition Am^ > is not satisfied due to the mass ordering of 
heavy Majorana neutrinos Eq. f fT8|) corresponding to light neutrino mass ordering. 

From Eqs. ([9]), (TT4]) and ( 1261) . the PMNS matrix in Eq. (l25l) can be expressed as 



PMNS 



v[ 22 v 21 -vi 23 v sl + vli 2 Vu v ^- y ^ +v ^ v[ 22 v 23 -v[ 23 v 33 + vti 2 v r3 

V V&^ss + V&&Vn + V[* 23 V 2 i YL^^lk V i 33 v 33 + vl hVl3 + F /* 3 y 23 y 

12 



P„(33) 



where is the (ij)-element of the mixing matrix V[, and Vij is the (r/')-element of U}JJ V 



given by 



v = ulu v 



( 



\ 



73 



1 

V3 



ce -(^l + f)_ se ^2-f) J_ 

^3 



73 

^Vl + f ) +ce '(V2-f ) 



)_ se '(-P2 + f ) j_ se '(Pl- j) +ee '(¥>2 + f) 

73 73 73 



\ 



(34) 



where s = sin# and c = COS6 1 . By recasting Eq. (133]) with the transformations e — >• e e iai , 
fi — y ji e l/?1 , t — y t e i/?2 and v 2 —y v 2 e l ^ ai ~ a2 \ we can rewrite the PMNS matrix as 



PMNS 



I C/ r ie-^ 2 U^ ai - aa -^ \U t3 \ J 



P' 

V 



(35) 



which corresponds to the standard parametrization as in PDG 20| . From the above equation, 
the neutrino mixing parameters can be displayed as 



sin 2 # 
sin 2 # 



\U, 



c2 



\u, 



12 



i - \u e3 \ 2 ' 



13 



\u, 12 



e3 1 j 



sin 2 9 2 3 

i-\u e3 \ 



(36) 



where «i = arg(C/ e i), «2 = arg([/ e2 ), a 3 = arg(C/ e3 ), fa = axg(C/ /i3 ) and /3 2 = arg(£/ T3 ). From 
the form of E/pmns given in Eq. (134j) . the solar, atmospheric and reactor neutrino mixing 
angles #12, #23 an d #13 can be approximated, up to order A 3 , as 

1 - 2 A cos <p e 3 + A 3 (cos (f) e 3 - 2 A cos <p e 2 + 2B cos $ 3 ) 



sin 2 # 



12 



3 - e - EX + A 2 # + £l\ 3 



2 „ 1 - sin 26 cos(^i2 - tt/3) -EX + A 2 T + A 3 6 
sin #23 = 



3 - e - EX + A 2 # + A 3 fi 



1 



sin# 13 = -^Ve + EX- A 2 ^ - X 3 Q , 
v3 

where 0f • = 0^ + 0^ and the parameters \l/ and Q are defined as 



(37) 



£ 

T 



1 + sin 29 cos <^i 2 , 

sin 2$ cos(<y?i2 — 7r/6) 



+ A 



v/3 cos 26* sin 2 ~~ cos 1 



2 sin 2$ cos((^i 2 — 7r/3)) 



7T 

cos 03 + sin 03 cos 2$ + 2 sin 26 cos 03 cos (^12 H — ) 

3 

\1> + 5 cos 0i (1 - 2 cos y?i2 sin 29) + sin 4>{ cos 20 



(38) 



13 



with 



e = = + B 



COS i 



4, 



1 + 2 sin 29 cos — tt/3)) - v / 3cos20sin^ 3 l . (39) 



In Eq. (1371) . if we turn off the higher dimensional operators in the Lagrangian, that is, if 
9 -> ±vr/4 and (p u -»■ 7r(0), and A -)• 0, the TBM angles sin 2 6>i 2 = l/3,sm 2 6> 23 = 1/2 and 
sin $i 3 = are restored, as expected. In the limit of 9 — > ir/A and (pu — > 7r (normal hierarchy 
of the neutrino masses), or — > —71/ 4 and <£> 12 — > (inverted hierarchy of the neutrino 
masses), the parameters behave as e — > 0, H — )■ 0, 'f — >■ —3/2, f2 — > 0, T — >■ — 3/2 — 3-B cos 0j 
and O — > 0. Then, the neutrino mixing angles can be simplified as 



1 2Acos0 3 ' . 2 , 1 . ( , t 1\ ,9 . « /A 2 , e + SA 



sin 2 e l2 « g-^, sin 2 ^3 ~ 2 + ( v cos ^ ~ 4 J A ' sin ^ 13 - V y + ~1T~ ' ^ 40) 

where A, B = 1 are used. Note here that in the limit of e — > and H — >■ the mixing param- 
eter sin^ goes to the value X/y/2. Leptonic CP violation at low energies can be detected 
through the neutrino oscillations which are sensitive to the Dirac CP-phase, but insensitive 



to the Majorana CP-phases in £/pmns [21] : the Jarlskog invariant Jcp = Im [6^1 6^2 ^2^1]' 
where U a j is an element of the PMNS matrix in Eq. (133 p . with a = e, /i, r corresponding to 
the lepton flavors and j = 1,2,3 corresponding to the light neutrino mass eigenstates. To 
see how both CP phases <y9 12 (coming from the neutrino sector) and 0f 2 3 (coming from the 
charged lepton sector) are correlated with low energy CP violation measurable through the 
neutrino oscillations, let us consider the leptonic CP violation parameter Jcp '■ 

cos 29 \\/3 f . 

Jcp = 7= — sin 29 sin 3 cos((^i 2 — n/6) 

6v3 9 

A 2 / 

cos 29 — B sin 29 sin <\> l x sin yj 12 
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sin 29 cos(</?i2 — 7r/6) sin0 3 + 2A sin 29 sin 2 cos(<^i2 + 7r/6) 
+ £>< 2 cos 2^ cos 0^ + sm 

20^ cos (/?i2 sin 

+ sinv9i 2 (cos03sin0^ + 3sin0^cos0i)j| + C(A 4 ) . (41) 

Note that sin 0g appears in the terms of the first order in A, and its value is crucial to satisfy 
the neutrino data for the solar mixing angle. In particular, for 9 — > ±7r/4 and <£> 12 — > 7r (or 0), 
we obtain Jqp — ±A/6 for sin 03 ~ ±1. 

In our numerical work, we use the five neutrino experimental data of Am 2 ol , Am 2 tm , 
#12, #13 and 6*23 at 3a level given in Eq. (JT]) as inputs. On this work, we are going to 
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FIG. 2: Left-plot displays the allowed regions of the solar mixing angle 9\2 versus the CP phase 
and right-plot displays the allowed regions of the atmospheric mixing angle #23 versus the CP 
phase The horizontal dotted lines in both plots correspond to the upper and lower bounds, 
respectively, satisfying the 3a experimental data in Eq. ([T]). 

study the case M 3 <C M 2 < Mi giving rise to a normal hierarchical mass ordering of light 
neutrinos, see Eq. f[3"2"j) . From Eqs. (I22|24|l28p we obtain 

m 2 = mo/(s.) - — ^ L ^jf L (42) 

which should be order of 0(0.01) eV, indicating Xp\y u \ 2 /M ~ C(10~ 15 ) GeV -1 . By using 
the relation m = Am 2 |^| 2 /167r 2 M with Eqs. ( I22ll42|) and the SM Higgs VEV v = 174 
GeV and, for example, by taking the scale M — 10 TeV and A*'' = 10 -8 , the values of the 
relevant parameters satisfying the five neutrino experimental data are taken as 



0.98 < k < 1.02 , 0.032 < \ yi/ \ < 0.045 , 0.0001 < y 1>2 < 0.47 
'1,2,3 



< <f>{ 2 3 < 2tt , < pi i2 < 2tt . (43) 



Without loss of generality, we take A = B = 1 appearing in the charged lepton mixing 
matrix V[, because they do not affect the leptonic mixing mixing parameters significantly. 
The left plot in Fig. [2] shows the allowed region of the solar mixing angle 812 versus the CP 
phase 0g. As can be seen the simplified equation of sin 2 9 12 in Eq. ( 140~1) . the mixing parameter 
#12 is dominantly controlled by the phased allowing 55° < <p\ < 100° and 260° < 0| < 300°. 
And the mixing parameter #23 is dominantly controlled by (f>\ as can be seen in Eq. (1401) . 
and the right plot in Fig. [2] displays the allowed region of the atmospheric mixing angle #23 
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FIG. 3: Allowed values for the reactor angle #13 in left-plot and Jcp in right-plot as a function of 
3 and #13, respectively. The horizontal and vertical dotted lines in the both figures correspond to 
the upper and lower bounds in 3<r of experimental data in Eq. ([1]). 

versus the CP phase cf)\. The left plot of Fig. |3] shows that the behavior of 6*13 as a function 
of S, where there is a lower bound 6* 13 > 3.5°. Because e > 0, depending on the sign of 
H, the second term in the squared-root of sin #13 in Eq. (HOT) can increase or decrease the 
value of #13 around the value Furthermore, because the value of H is bounded by 

as can be seen in Eq. ( 138]) . we expect that there is a lower bound on the possible value of 
6*13. The parameter H depends mainly on y\ and 1/2, defined in Eq. ([8]), which represent the 
effects of the dimension-5 operators. Thus, in our scenario the lower bound on the mixing 
angle 6*13 is strongly dependent on the cutoff scale A, the A4 symmetry breaking scale v x 
and the relevant couplings \y% a \, through yi and y 2 . Since neutrino oscillation experiments 
are sensitive to the Dirac CP phase Sep, the Jarlskog invariant of the leptonic sector given 
in Eq. (l4Tj) would be a signal of CP violation. The right plot in Fig. [3] shows that allowed 
values for Jcp as a function of #13, showing \Jcp\ ~ 0.01 — 0.04 due to the sizable #13. 
This can be tested in the future experiments such as the upcoming long baseline neutrino 
oscillation ones. In the plots Fig. [2] and Fig. |3] the horizontal and vertical dotted lines in 
the both figures correspond to the upper and lower bounds in 3a of experimental data in 
Eq. ©. 
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IV. DARK MATTER AND LEPTON FLAVOR VIOLATION 



In this section, we study a fermionic DM analysis, that is, right-handed neutrinos. Espe- 
cially, we are interested in a hierarchical case of heavy Majorana neutrino, M 3 <C M 2 < Mi 
for £ = 0, giving a hierarchical normal mass ordering of light neutrinos, as shown in section- 
Ill. 

The existence of the flavor neutrino mixing implies that the individual lepton charges, 



L a ,a = e,fj,,T are not conserved [15| and processes like £ a — > £37 should take place. Ex- 
perimental discovery of lepton rare decay processes £ a —> £37 is one of smoking gun signals 
of physics beyond the SM; thus several experiments have been developed to detect LFV 



processes. The present experimental upper bounds are given at 90% C.L. [22J, |23j as 

Br(/i ->■ erf) < 2.4 x 10~ 12 , Br(r ->■ ^7) < 4.4 x 10~ 8 , 

Br(r -> ery) < 3.3 x 1(T 8 . (44) 

Assuming the lightest particle of Ni to be DM, we consider annihilation of Ni through 
Yukawa interaction Eq. ( fl~5|) in the early Universe. The DM mass Mii g ht es t — M% is con- 
strained by the LFV processes, if the size of \y u \ which is the overall scale of Y v is fixed 
through the DM relic density. For the given Yukawa interaction of the fermion singlet with 
SM particles, Y^i^Ni, its annihilation rate into the latter, as shown in Fig. HI and its relic 
density can be calculated and are related to each other by the thermal dynamics of the 



Universe within the standard big-bang cosmology [24j. Then in the WMAP analysis we find 
the thermally averaged cross section (av) for the annihilation of two iV 3 's from Fig. H]in the 
limit of the vanishing final state lepton masses: 

(av) c '^ |V f ;J + 2r2) , r = M 3 2 /K ± + Mil (45) 

where we assume m^o = m rj ±. The thermally averaged cross section Eq. (j4"5"]) does not contain 
s-wave contribution as a consequence of massless limit of the final state particles. Assuming 
this to be the dominant contribution to the DM relic density of the Universe, we need that 
the pair annihilation cross section must be of the order of 

(a eS v r ) ^ 8.35 x 10~ W ^L ( ' GeV~ 2 , (46) 



fg- V 0.11 

which is the thermally averaged annihilation cross section of a pair of fermion particles into 
SM particles and the relative speed of the fermion pair in their center-of-mass frame, where 
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FIG. 4: Annihilation diagrams of n for the cross section av, where and v run over three families. 



the decoupling temperature Tf of particles, xj ~ In Q.038Mp ^nja off |iy| ) ^ ^ e curren t Particle 



Data Group value for the DM density, Q d h 2 = 0.1123 ± 0.0035 [16|, in which Q d is the 
cosmological parameter associated with DM, h is the normalized Hubble constant, and 
g* is the number of relativistic degrees of freedom at the DM chemical coupling. From 
Eqs. ( )45|46|) . one can obtain for m v ± = x v GeV 

- - 
2 / (av) \ 4 / m, n ± \ 2 



\yj ~ 4.3 x 10- 2 VxT i— i , J> - , (47) 

V " V3 x 10- 9 GeV"V \x v GeVJ V ; 

where r ~ 1/2 is used. From Eq. (|4"7]) we see that the coupling depends on the scale of dark 
matter. For the weak coupling to be ensured, one needs \y v \ < 4ir 3 which, in turn, indicates 
the upper bound on the DM mass must be roughly 

M 3 < 83 TeV . (48) 

Now let us consider LFV, especially r — > /ry for a lower bound of DM and \x — > e'-f for a 
connection to the neutrino parameter # 13 . Fig. depicts that one-loop diagrams to the one 
for neutrino masses contribute to the lepton flavor violating processes like £ a —> (a, (3 = 



e, t), whose branching ratio is estimated as [26] 

Q 

Br(i a £ 01 ) = 647r(G ,"^ 2±)2 i^l 2 Br(4 -> ^P^i/ Q ) (49) 
where a e ~ 1/137 and is the Fermi constant, and B a p is given by 

3 

B aP = Y,y ai Y; i F 2 {x l ) , (50) 



\y v \ < 1 case has been considered in Ref. 
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FIG. 5: One-loop diagrams to the one for neutrino masses contribute to the lepton flavor violating 
processes like l a — > £97 (a, f3 = e,fj,,r). 



in which ^(xj) is given by 

_ 1 - 6xi + 3x 2 + 1x\ - QxjhiXi 

F2[Xt) ~ 6(1 -x t y (51) 

withF2(l) = l/12andxj = M 2 /m 2 ± . Note here that B a p does not depend on the phases ipi^ 
associated with heavy Majorana neutrino mass matrix. From Eq. ( 1491) the LFV branching 
ratios can be simplified as 



\B, 



12 

Br(4 ~ 8.0 x 10 5 r a/3 ^^GeV 4 , (52) 

m r,± 

where r Me = 1.0, r Te = 0.1784 and r TfJi = 0.1736, which indicates that the LFV branching 
ratios depends on m^i and on the neutrino parameters in a flavor dependent manner. 

In the following we consider the case, M 3 ~ m v ± <C M 2 < Mi, as shown in section-Ill. 
Taking the case M 1 ~ 2M, M 2 ~ M and M 3 ~ 0.01M into account, the function F 2 (xi) can 
have the values 1/12 (x 3 = 1), 3.3xl0~ 5 (x2 = 10 4 ) and 8.3 x 10" 6 ( x l =4xl0 4 ), for N U N 2 
and N 3 , respectively, which indicates only the lightest heavy Majorana neutrino among the 
heavy neutrinos can contribute the branching ratio of LFV. Then, the expressions of |-B a /?| 2 
relevant for r — > ^7, r — > e'-f and /i — > e7, respectively, are approximately given as 

\B T „\ 2 = \y u \ 4 F 2 (x 3 )Q Tfl , \B Te \ 2 ~ |£? Me | 2 = \y u \ 4 F 2 (x 3 )Q^ e , (53) 

where the quantities Q T ^ and Q^ e are simply expressed up to order of A 2 as 

2 



1 2/1 cos pi +2/2 cos p 2 A , . t e s A 

Qt^ - j 2 ^/f ^ 2 sm ^ 2 ~ ^3) - yi sm (pi - ^3)) - -j 

Q, e - (yi 7/ 2)2 +^sin 2 ^ 1 - p2 



■ 



12 3 V 2 

A A 2 

(j/ 2 sin(p 2 - 4&) - Vi sin(pi - 3 )) + — . (54) 
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FIG. 6: The Q^ e < x 1.26 x 10 10 sensitivity to #13. Here the vertical dotted lines represent 
the upper and lower bounds in 3<r of Eq. ([!]) 



From the constraint of the branching ratio of r — > /ry in Eq. (1441) we can obtain a lower 
bound of DM mass by using Eqs. A47I52|54[) : 

66 GeV < M 3 < 83 TeV , (55) 

where the upper bound is from Eq. (I48p . Using Eq. (jHJ), the branching ratio Eq. (152]) can 

be expressed in terms of Q^ e , Q Tfl and x v as 

Que < x\ x 1.26 x 1(T 10 , Q ni <x 2 v x 1.33 x 10" 5 , (56) 

indicating Q^ e is much more stringent than Q Tii in sensitivity. In Fig. [6j the allowed region 
between 8 13 and is described. It tells us that a rather large #13 > 9° is predicted for 12.6 
TeV < M3 < 14.6 TeV. For M3 < 6.3 TeV, on the other hand, a rather smaller 6*13 < 8° is 
predicted. It might be able to detected by direct or indirect DM searches, though the mass 
sale is rather high. 

Finally we mention the detectability at the LHC. We find that an inert charged Higgs, 
which decays into a muon and tauon (or a anti-muon and a anti-tauon) universally with a 
large missing energy in the limit of leading order, where the missing energy is carried away 
by our DM. It suggests that the clean signal will be obtained at the LHC. 
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V. CONCLUSION 



We have analyzed the lepton masses and the mixings based on A4 flavor symmetry. 
Especially, we obtained a lower bound of 6* 13 > 3.5° in hierarchical normal mass ordering. 
And, we found that the large #13 can be correlated with the Majorana type of DM, and its 
mass be (9(1-10) TeV, as a result of analyzing the WMAP and lepton flavor violation such 
as the \i — > ej experiment. Moreover, we could predict that a rather lager # 13 ; 9° < #13 for 
12.6 TeV < M 3 < 14.6 TeV, and a rather smaller 9 13 ; 9 13 < 8°, for M 3 < 6.3 TeV. It implies 
that it might be able to detetected by direct or indirect DM searches as well as LHC at a 
rather high energy scale. 



Appendix A: Higgs Potential and vacuum alignment 

Since it is nontrivial to ensure that the different vacuum alignments of (tp°) = (v, v, v) and 
(x) = i v xi 0' 0) i n Eq. ([6]) are preserved, we shall briefly discuss these vacuum alignments. 
There is a generic way to prohibit the problematic interaction terms by physically separating 
the fields x an d Here we solve the vacuum alignment problem by extending the 

model with a spacial extra dimension y 6j. We assume that each field lives on the 4D 
brane either at y = or at y = L, as shown in Fig. [3 The heavy neutrino masses arise 
from local operators at y — 0, while the charged fermion masses and the neutrino Yukawa 
interactions are realized by non-local effects involving both branes. A detailed explanation 
of this possibility is beyond the scope of this paper. 



Qr Ql 

Ir L l 

X $ 

Nr V 

y l 

FIG. 7: Fifth dimension and locations of scalar and fermion fields. 

Then, the most general renormalizable scalar potential of $,77 and \i invariant under 
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SU(2) L x U(1) Y x A 4 x Z 2 , is given by 

Vy= L = + Af ($t$) 1 ($t$) 1 + A|($t$) 1 ,($t$) 1 „ + Af ($t$) 3s ($t$) 3s 

+ Af ($t$) 3a ($t$) 3a + A * ($ t $)3s($ t $)3a + + A , ( ^ )2 

+ Af^^^irt) + \p(&r])(rf$) + A* ,7 ($ t ^)($ t ?7) + Xf^tftytf®) > ( A1 ) 
1^=0 = + £ x (xxx)i + A?(xx)i(xx)i + >$(xx)i'(xx)i» + A3 (xx)s(xx)3 

+ ^{xx)i{xx)sx + -r(xx)i'(xx)3X + ^(xx)i"(xx)sx + ^(xx)3(xx)3X,(A2) 

where //,,, \x x and £ x are of the mass dimension 1, while Af 5 , X , A x 3 , Af^ 3 and «i,...,4 
are all dimensionless. Notice also here that there are no five dimensional terms in V v =l and 



K 2 = K 3 = 0. Let us set tp® 



ip Ri + i (pn (i = 1, 2, 3) and rp = rp R + i rf R . Then masses 



of $j, 77, and Xi are written as follows: 



m\ = m 2 (Re(?7)) 



1 



m, 



m 



2 (Mv)) = \4 + r 2 ( X i V + A 2 ^-2A^) 



2 1 2 2 \ ^T) 



2' " 2 

/ 2(Af + Af) 2Af - Af + 4Af 2Af - Af + 4Af \ 
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Re Mi 



2Af 



Af + 4Af 
1* 



2(Af + Af) 2Af-Af + 4Af 



/ 
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Im(i/))i 



-2(Af + A*) 
(Af + A*) 
(Af + A*) 



-2A 3 A 3 



(Af + Af) (Af + Af) N 
-2(Af + Af) (Af + Af) 
(Af + A*) -2(Af + A*) J 



Af 
Af 



-2Af 
Af 



A 3 

Af 
-2Af 



I 2Af - Af + 4Af 2Af - Af + 4Af 2(Af + Af ) J 



(A3) 
(A4) 
(A5) 

(A6) 



(A7) 



(A8) 



rn 



[RcM-ImMli 



~ A 5 A 5 



A 5 







_A 5 A 5 



_A 5 




/ 



m. 



(A9) 



. Wi = rfiay ^2(AJ + A*) , A? - ^A* + 2A* , A? - ^ + 2A^ , (A10) 
where we assume Af^ is real to forbid the mixing term between Re(r/) and Im(r/), and also 
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X is real. Here we use stable conditions; /i| = — (6Af + 8Af )v 2 and = — 2(A* + A^)^ 
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